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An algorithm of Monte Carlo code for the production of fast hadron 
I/") ■ leader in the proton-(anti)proton interaction by QCD process is dis- 

cussed. 
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1 Introduction 

At present one of a few most interesting questions is theoretical explanation of 
the experimental data for hevi-ion collisions. Therefore any alternative expla- 
nations of phenomena which could be interpreted as a signature of new physics 
such as quark-gluon plasma deserve careful elaboration. In the Ref. we place 
for consideration a new hadronization process based on the cascade decay of 
rotating strings into hadron states. As it had marked in Ref.|l|, long ago it 
have been already known fact that linearity of Regge trajectory could be ex- 
plained by supposition that the hadron is a rigid straight rotating string |2|]. If 
the Regge trajectories are unlimited, it gives a chance to consider the rotating 
strings having large masses too. There is no a complete theory for decay of 
rotating string, and so we have been constructing a Monte Carlo parametriza- 
tion for this process [Q]. In our cascade string breaking model the process of 
breaking of mother string into two daugher strings is repeated. The string 
rotation leads to correlation between the transverse momentum of daughter 
string and the mass of mother string. Therefore the average transverse mo- 
mentum of secondary hadrons, produced by decays of rotating strings, grows 
with the mass of primary string (see Fig.5 in |l| and Fig. 10 in ||), and so 
the fluctuations of the masses of primary strings (in the model for proton- 
antiproton collisions) lead to explanation of experimental Pt distribution of 
secondaries up to pt = 4 Gev at 1800 GeV for proton-antiproton interactions 
(see Fig. 13 in Ref.|| where there are preliminary results). 

For heavy-ion reactions there is widely known theoretical prediction that 
hard parton-parton interaction and formation of quark-gluon plasma are two 
sources of secondary particles with high transverse momenta. We hope that 
our preliminary calculations shown that the string rotation is also the source of 
high transverse momenta, and so the string rotation can lead to effect which is 
the very same as vapourization of drops of quark-gluon plasma. Therefore this 
result can't be ignored in the realistic calculations of high energy reactions. 

The starting point for nucleus-nucleus interactions is nucleon-nucleon in- 
teraction, and so we should be constructed proton- (anti)proton Monte Carlo 
model for high energy interaction where, according to cosmic ray experimen- 
tal data[|], ||, p, for interaction at energy E Lab > 10 16 eV there is a group of 
very fast hadron leaders, which have not till now the satisfactory theoretical 



explanation^ |5], |6|. 

For the modelling calculations of production of the fast hadron leaders 
at high-energy proton- (anti)proton interactions we suppose to use a parton 
model idea that (anti)proton is a collection of quasi-free partons which share 
its momentum, and that the partons carry negligible transverse momentum. 
In many models of multiparticle production (see, for example, J7|, |j| ) in the 
first approximation^ the parton structure of proton is a color-triplet quark 
and an antitriplet diquark, but in the standard parton model it is many 
(anti)quarks and gluons^, ||. Therefore in our calculations we should be 
taken into account that each incoming hadron is composite object, consisting 
of constituent (anti)(di) quarks which have possibility to dissociate into quasi- 
free (anti)(di) quarks and gluons (constituent partons). This process can lead 
to formation of excited partons by interaction of partons of target and pro- 
jectile (for example, a gluon of target can be "absorbed" by (anti)(di) quark 
of projectile). Invariant mass of excited (anti)(di) quark can be choosed as the 
characteristic of excitation. If this excitation is sufficient to starting of QCD 
cascade, it should be calculated. If this excitation is under threshold of QCD 
cascade, we can generate "decay" of excited parton into non-excited parton 
and hadron leader (see Fig. la). Therefore for the mass of excited parton there 
is a threshold which is parameter of modelling and which is a minimal in- 
variant mass sufficient to start of QCD cascade, and simultaneously it may 
be a maximal invariant mass of excited parton sufficient to start of hadron 
leader production. However, it should be marked, there is not any restriction 
to production of the fast hadron leader by excited parton with invariant mass 
greater of threshold sufficient to start of QCD cascade which has priority at 
high invariant mass of excited parton. Of course, fast hadron leader can be 
produced also after QCD cascade when invariant mass of cascading parton 
becomes less than the value of threshold. 

If (constituent) (di)quark and anti(di)quark can't be decayed, they stretch 
the primary string which is decayed into secondary hadrons according to the 
algorithm described in M. 



1 In the frame of our model a qualitative explanation for cosmic ray experimental data 
about alignment of hadron leadersW, [5| |(| can be given by additional physics connected with 
the fast hadron leader production described here and with cascade decay of rotating strings 
into hadron statesllj. However, it requires detailed description which overstep the limits of 
the theme of presented paper. Therefore the detailes will be given elsewhere. 

2 Only valence quarks are considered. 



2 Formation of excited partons 

An algorithm of Monte Carlo (MC) code of formation of excited partons is 
described in this Section. In our model (as in Ref.fllTfl) constituent quark of 
proton can dissociate into quasi-free quark and gluon which has the part y 
of four-momentum of constituent quark. The probability of dissociation into 
interval dy is formally determined by perturbative formula (see Appendix 1 
where the variable z — 1 — y) 

1 + (1 _y) 2 

dW = a dy at y min < y < 1 , (1) 

y 

and integral probability is determined by the formula 

W = a i ^- dy at y min < y < 1 , (2) 



ymin y 



* 1 + (1 - yf ,_! ( . 

a = [ dy } , (3) 



where 



'2/i y 

ymin = |/i at 2/i > 2/o = (£o/xE) , (4) 

y m in = 2/0 at 2/1 < 2/o < 1 , (5) 

2/mm = 1 at 1 < 2/0 , (6) 

Eq and ?/i are free parameters of modelling, x is part of four-momentum which is 
carried by constituent (anti)(di) quark before dissociation of it, E is the energy 
of (anti)proton in c.m.system of collising particles. For example, the values 
can be choosed as 

e ~ m u , y 1 < 0.01 , (7) 

where m u = 0.34GeV is the mass of constituent w-quark. As one can see from 
eqs.(2)-(6), the probability of dissociation grows with energy xE of constituent 
quark. 

The parts x and %' of four-momentum are determined by parton distribu- 
tion function for the projectile and target. If (according to eq.(2)) there is the 
process of dissociation of constituent (anti)(di) quark into (anti)(di) quark and 



gluon, the gluon part y of four-momentum of constituent quark is generated 
by the formula 

• ^C-y) 2 dy i f 1 i±^U . (8) 



Vmin y J Vmin y 

where r is uniformly distributed random number from the range (0,1). By 
analogy with eq.(8) the part y' for gluon from the target can be generated. In 
this case, in eq.(4) instead y the value y' = (eq/x'E) should be calculated. 
Let us now to consider an example for the formation of excited parton q'* 
by interaction 

g + Y c - 9* , (9) 

i.e. a gluon of target is "absorbed" by antiquark of projectile at proton- 

antiproton collision?]. 

Before simulation of the subprocess (9), the probability of this process 

should be determined. Let p is density of current of protons to antiprotons 

and d 3 p is given at fixed flowers of constituents q c and q' c from the intervals 

dx, dx', dy provided that constituent antiquark q' c is not dissociated; and let 

o is cross-section of inelastic non-diffractive pp collisions , at is cross-section 

of the subprocess (9) , then the probability of the subprocess (9) is given by 

the formula 

d 3 pat 



a- 



W* = —J^ = -T • (10) 



q p d 3 a a 

In the eq.(10) it was taken into account that 

d 3 p d 3 cr 

p a 

According to dimension of cross-section we can choose 



(11) 



at oc —, (12) 



M*f 



for the invariant mass 



Mlf = (x 2 y 2 + x' 2 )M 2 + 2xyx'(E 2 + P 2 ) (13) 



3 It should be marked, the subprocess q c + q' c —> (q + g) + (q' + g') — > (q + q') + (g + g') 
9i + 32 can't be a source to production of hadron leader. 



which is under threshold sufficient to start of QCD cascade, or 

<% « "4^ (14) 



Q 

for the invariant mass sufficient to start of QCD cascade^. In eqs.(12)-(14) Mt 
is invariant mass of excited antiquark q'* from (9) and P is the momentum 
of (anti)proton in cm. system of collising particles, a s {M^?) is the strong cou- 
pling constant. However, the probability value should be at least restricted. 
Therefore for small invariant mass which is under threshold sufficient to start 
of QCD cascade, we can choose an approximation of probability for the process 
(9) as 

w i « JET^ ' (15) 

where (3 is free parameter of modelling. 

The eqs.(2)-(8), (10), (13)-(15) are sufficient for the modelling of starting- 
point for QCD cascade or production of fast hadron leader. 



3 Fast hadron leaders 

Short distance QCD interaction leads to formation of excited (anti)(di) quark 
state with short life time which is less than time sufficient for the forma- 
tion of long distance string. Therefore the process of "decay" of excited 
(anti)(di) quark with invariant mass under threshold sufficient to start of QCD 
cascade , i.e. the process of fast hadron leader production (see Fig. la) is before 
formation of string. 

In Fig. la it is shown a "decay" of excited quark q*. Therefore this diagram 
is related with QCD diagram in Fig. lb where a quark q is carried the part of 
four-momentum z of excited quark q* distributed (see Appendix 1) by the law 

1 + z 2 
f(z) oc -^— . (16) 

Therefore in Fig. la we can choose for the "decay" of excited (anti)(di) quark q* 
the same variable z and law (16) which has peak in the point z—1 , i.e. hadron 

4 In ref.fll0| there are examples of cross-sections for the more complicated subprocesses. 



h is carried the part of four-momentum z of excited (anti)(di) quark q* near 
by 1 , i.e. this hadron h is fast. In the Fig. la slow quark q stretchs the string 
which is decayed into secondary hadrons according to the algorithm described 
in U||]. Therefore in Fig. la the hadron h is fast leader. 

In the frame of our model, each proton-(anti)proton interaction leads to 
production at least of one fast hadron leader which carrys large part of energy 
of primary particle (see eq.(16)). This process is repeated at second inter- 
action (in the atmosphere), and so, in our model , it can be understood an 
experimental phenomena that hadron leader has small dissipation of energy at 
interactions in the atmosphere 0. 

Appendix 1 

Let us remind of some results of QCDpT| which are important under con- 



struction of our model. We take a look at current of virtual photons 7* to 

quarks q (Fig. 2) at fixed cm. energy a/s. Let this current of virtual photons 
have uniform distribution in variable 

- Q2 , (17) 



8 + Q 2 

where Q 2 is squared four-momentum of photon taken with opposite sign. Dif- 
ferential cross-section of the process given by Fig. 2 is 

d 2 a = o\ "f2(PTi z ) &Vt dz at — t <C s , (18) 

where 

a s 1 4 1 + z 2 , s 

72 = it -2 o 1 , 19) 



o\ is cross-section of the block 1 in Fig.2 , a s is the strong coupling constant, pt 
is gluon transverse momentum given relative to momentum vector of primary 
quark in c.m.s. , t is the squared difference between the momenta of primary 
quark and gluon (see Fig.2). It is important , at —t<^s an intermediate 
quark carrys the part z of four-momentum of primary quark in c.m.s. 
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Figure Captions 



Fig.l a. The "decay" of excited parton q* into non-excited parton q and 
hadron leader h . b. QCD diagram for the decay of excited quark q* into 
quark q and gluon g . 

Fig. 2 One in a few diagrams of the process 7* q — > q g . 



This figure "FIGURE1.GIF" is available in "GIF" format from: 
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